Investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment.

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

Epigenomic mapping reveals distinct B cell acute lymphoblastic leukemia chromatin architectures and regulators.

B cell lineage acute lymphoblastic leukemia (B-ALL) is composed of diverse molecular subtypes, and while transcriptional and DNA methylation profiling has been extensively examined, the chromatin landscape is not well characterized for many subtypes. We therefore mapped chromatin accessibility using ATAC-seq in primary B-ALL cells from 156 patients spanning ten molecular subtypes and present this dataset as a resource. Differential chromatin accessibility and transcription factor (TF) footprint profiling were employed and identified B-ALL cell of origin, TF-target gene interactions enriched in B-ALL, and key TFs associated with accessible chromatin sites preferentially active in B-ALL. We further identified over 20% of accessible chromatin sites exhibiting strong subtype enrichment and candidate TFs that maintain subtype-specific chromatin architectures. Over 9,000 genetic variants were uncovered, contributing to variability in chromatin accessibility among patient samples. Our data suggest that distinct chromatin architectures are driven by diverse TFs and inherited genetic variants that promote unique gene-regulatory networks.

Epigenomic mapping in B-cell acute lymphoblastic leukemia identifies transcriptional regulators and noncoding variants promoting distinct chromatin architectures.

B-cell lineage acute lymphoblastic leukemia (B-ALL) is comprised of diverse molecular subtypes and while transcriptional and DNA methylation profiling of B-ALL subtypes has been extensively examined, the accompanying chromatin landscape is not well characterized for many subtypes. We therefore mapped chromatin accessibility using ATAC-seq for 10 B-ALL molecular subtypes in primary ALL cells from 154 patients. Comparisons with B-cell progenitors identified candidate B-ALL cell-of-origin and AP-1-associated cis-regulatory rewiring in B-ALL. Cis-regulatory rewiring promoted B-ALL-specific gene regulatory networks impacting oncogenic signaling pathways that perturb normal B-cell development. We also identified that over 20% of B-ALL accessible chromatin sites exhibit strong subtype enrichment, with transcription factor (TF) footprint profiling identifying candidate TFs that maintain subtype-specific chromatin architectures. Over 9000 inherited genetic variants were further uncovered that contribute to variability in chromatin accessibility among individual patient samples. Overall, our data suggest that distinct chromatin architectures are driven by diverse TFs and inherited genetic variants which promote unique gene regulatory networks that contribute to transcriptional differences among B-ALL subtypes.

Functional investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment.

Although acute lymphoblastic leukemia (ALL) is the most common childhood cancer, there is limited understanding of the contribution of inherited genetic variation on inter-individual differences in chemotherapy response. Defining genetic factors impacting therapy failure can help better predict response and identify drug resistance mechanisms. We therefore mapped inherited noncoding variants associated with chemotherapeutic drug resistance and/or treatment outcome to ALL cis-regulatory elements and investigated their gene regulatory potential and genomic connectivity using massively parallel reporter assays and promoter capture Hi-C, respectively. We identified 53 variants with reproducible allele-specific effects on transcription and high-confidence gene targets. Subsequent functional interrogation of the top variant (rs1247117) determined that it disrupted a PU.1 consensus motif and PU.1 binding affinity. Importantly, deletion of the genomic interval containing rs1247117 sensitized ALL cells to vincristine. Together, these data demonstrate that noncoding regulatory variation associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to chemotherapeutic agents in ALL.

Human SMARCA5 is continuously required to maintain nucleosome spacing.

Genetic models suggested that SMARCA5 was required for DNA-templated events including transcription, DNA replication, and DNA repair. We engineered a degron tag into the endogenous alleles of SMARCA5, a catalytic component of the imitation switch complexes in three different human cell lines to define the effects of rapid degradation of this key regulator. Degradation of SMARCA5 was associated with a rapid increase in global nucleosome repeat length, which may allow greater chromatin compaction. However, there were few changes in nascent transcription within the first 6 h of degradation. Nevertheless, we demonstrated a requirement for SMARCA5 to control nucleosome repeat length at G1/S and during the S phase. SMARCA5 co-localized with CTCF and H2A.Z, and we found a rapid loss of CTCF DNA binding and disruption of nucleosomal phasing around CTCF binding sites. This spatiotemporal analysis indicates that SMARCA5 is continuously required for maintaining nucleosomal spacing.

Mutual antagonism between glucocorticoid and canonical Wnt signaling pathways in B-cell acute lymphoblastic leukemia.

Glucocorticoids (GCs; i.e., steroids) are important chemotherapeutic agents in the treatment of B-cell precursor acute lymphoblastic leukemia (B-ALL) and de novo GC resistance predicts relapse and poor clinical outcome in patients. Glucocorticoids induce B-ALL cell apoptosis through activation of glucocorticoid receptor (GR), a ligand-induced nuclear receptor transcription factor (TF). We previously identified disruptions to glucocorticoid receptor (GR)-bound cis -regulatory elements controlling TLE1 expression in GC-resistant primary B-ALL cells from patients. TLE1 is a GC-response gene up-regulated by steroids and functions as a canonical Wnt signaling repressor. To better understand the mechanistic relationship between GC signaling and canonical Wnt signaling, we performed diverse functional analyses that identified extensive crosstalk and mutual antagonism between these two signaling pathways in B-ALL. We determined that crosstalk and antagonism was driven by the binding of GR and the canonical Wnt signaling TFs LEF1 and TCF7L2 to overlapping sets of cis -regulatory elements associated with genes impacting cell death and cell proliferation, and was further accompanied by overlapping and opposing transcriptional programs. Our data additionally suggest that cis -regulatory disruptions at TLE1 are linked to GC resistance through a dampening of the GC response and GC-mediated apoptosis via enhanced canonical Wnt signaling. As a result of the extensive genomic and gene regulatory connectivity between these two signaling pathways, our data supports the importance of canonical Wnt signaling in mediating GC resistance in B-ALL.

Epigenetic activation of the FLT3 gene by ZNF384 fusion confers a therapeutic susceptibility in acute lymphoblastic leukemia.

FLT3 is an attractive therapeutic target in acute lymphoblastic leukemia (ALL) but the mechanism for its activation in this cancer is incompletely understood. Profiling global gene expression in large ALL cohorts, we identify over-expression of FLT3 in ZNF384-rearranged ALL, consistently across cases harboring different fusion partners with ZNF384. Mechanistically, we discover an intergenic enhancer element at the FLT3 locus that is exclusively activated in ZNF384-rearranged ALL, with the enhancer-promoter looping directly mediated by the fusion protein. There is also a global enrichment of active enhancers within ZNF384 binding sites across the genome in ZNF384-rearranged ALL cells. Downregulation of ZNF384 blunts FLT3 activation and decreases ALL cell sensitivity to FLT3 inhibitor gilteritinib in vitro. In patient-derived xenograft models of ZNF384-rearranged ALL, gilteritinib exhibits significant anti-leukemia efficacy as a monotherapy in vivo. Collectively, our results provide insights into FLT3 regulation in ALL and point to potential genomics-guided targeted therapy for this patient population.

Sustained Minimal Residual Disease Negativity in Multiple Myeloma is Associated with Stool Butyrate and Healthier Plant-Based Diets.

PURPOSE: Sustained minimal residual disease (MRD) negativity is associated with long-term survival in multiple myeloma. The gut microbiome is affected by diet, and in turn can modulate host immunity, for example through production of short-chain fatty acids including butyrate. We hypothesized that dietary factors affect the microbiome (abundance of butyrate-producing bacteria or stool butyrate concentration) and may be associated with multiple myeloma outcomes. EXPERIMENTAL DESIGN: We examined the relationship of dietary factors (via a food frequency questionnaire), stool metabolites (via gas chromatography-mass spectrometry), and the stool microbiome (via 16S sequencing - α-diversity and relative abundance of butyrate-producing bacteria) with sustained MRD negativity (via flow cytometry at two timepoints 1 year apart) in myeloma patients on lenalidomide maintenance. The Healthy Eating Index 2015 score and flavonoid nutrient values were calculated from the food frequency questionnaire. The Wilcoxon rank sum test was used to evaluate associations with two-sided P < 0.05 considered significant. RESULTS: At 3 months, higher stool butyrate concentration (P = 0.037), butyrate producers (P = 0.025), and α-diversity (P = 0.0035) were associated with sustained MRD negativity. Healthier dietary proteins, (from seafood and plants), correlated with butyrate at 3 months (P = 0.009) and sustained MRD negativity (P = 0.05). Consumption of dietary flavonoids, plant nutrients with antioxidant effects, correlated with stool butyrate concentration (anthocyanidins P = 0.01, flavones P = 0.01, and flavanols P = 0.02). CONCLUSIONS: This is the first study to demonstrate an association between a plant-based dietary pattern, stool butyrate production, and sustained MRD negativity in multiple myeloma, providing rationale to evaluate a prospective dietary intervention.

Epigenomic profiling of glucocorticoid responses identifies cis-regulatory disruptions impacting steroid resistance in childhood acute lymphoblastic leukemia.

Glucocorticoids (GCs) are a mainstay of contemporary, multidrug chemotherapy in the treatment of childhood acute lymphoblastic leukemia (ALL), and resistance to GCs remains a major clinical concern. Resistance to GCs is predictive of ALL relapse and poor clinical outcome, and therefore represents a major hurdle limiting further improvements in survival rates. While advances have been made in identifying genes implicated in GC resistance, there remains an insufficient understanding of the impact of cis-regulatory disruptions in resistance. To address this, we mapped the gene regulatory response to GCs in two ALL cell lines using functional genomics and high-throughput reporter assays and identified thousands of GC-responsive changes to chromatin state, including the formation of over 250 GC-responsive super-enhancers and a depletion of AP-1 bound cis-regulatory elements implicated in cell proliferation and anti-apoptotic processes. By integrating our GC response maps with genetic and epigenetic datasets in primary ALL cells from patients, we further uncovered cis-regulatory disruptions at GC-responsive genes that impact GC resistance in childhood ALL. Overall, these data indicate that GCs initiate pervasive effects on the leukemia epigenome, and that alterations to the GC gene regulatory network contribute to GC resistance.

Amino acid stress response genes promote L-asparaginase resistance in pediatric acute lymphoblastic leukemia.

Understanding the genomic and epigenetic mechanisms of drug resistance in pediatric acute lymphoblastic leukemia (ALL) is critical for further improvements in treatment outcomes. The role of transcriptomic response in conferring resistance to l-asparaginase (LASP) is poorly understood beyond asparagine synthetase (ASNS). We defined reproducible LASP response genes in LASP-resistant and LASP-sensitive ALL cell lines as well as primary leukemia samples from newly diagnosed patients. Defining target genes of the amino acid stress response-related transcription factor activating transcription factor 4 (ATF4) in ALL cell lines using chromatin immunoprecipitation sequencing (ChIP-seq) revealed 45% of genes that changed expression after LASP treatment were direct targets of the ATF4 transcription factor, and 34% of these genes harbored LASP-responsive ATF4 promoter binding events. SLC7A11 was found to be a response gene in cell lines and patient samples as well as a direct target of ATF4. SLC7A11 was also one of only 2.4% of LASP response genes with basal level gene expression that also correlated with LASP ex vivo resistance in primary leukemia cells. Experiments using chemical inhibition of SLC7A11 with sulfasalazine, gene overexpression, and partial gene knockout recapitulated LASP resistance or sensitivity in ALL cell lines. These findings show the importance of assessing changes in gene expression following treatment with an antileukemic agent for its association with drug resistance and highlight that many response genes may not differ in their basal expression in drug-resistant leukemia cells.

Dual detection of chromatin accessibility and DNA methylation using ATAC-Me.

The epigenome is multidimensional, with individual molecular components operating on different levels to control transcriptional output. Techniques that combine measurements of these features can reveal their precise correspondence in genomic space, or temporal connectivity, to better understand how they jointly regulate genes. ATAC-Me is an integrated method to probe DNA methylation and chromatin accessibility from a single DNA fragment library. Intact nuclei undergo Tn5 transposition to isolate DNA fragments within nucleosome-free regions. Isolated fragments are exposed to sodium bisulfite before library amplification and sequencing. A typical ATAC-Me experiment detects ~60,000-75,000 peak regions (P < 0.05), covering ~3-4 million CpG sites with at least 5× coverage. These sites display a range of methylation values depending on the cellular and genomic context. The approach is well suited for time course studies that aim to capture chromatin and DNA methylation dynamics in tandem during cellular differentiation. The protocol is completed in 2 d with standard molecular biology equipment and expertise. Analysis of resulting data uses publicly available software requiring basic bioinformatics skills to interpret results.

Lysine acetylation regulates the RNA binding, subcellular localization and inclusion formation of FUS.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease characterized by the preferential death of motor neurons. Approximately 10% of ALS cases are familial and 90% are sporadic. Fused in sarcoma (FUS) is a ubiquitously expressed RNA-binding protein implicated in familial ALS and frontotemporal dementia (FTD). The physiological function and pathological mechanism of FUS are not well understood, particularly whether post-translational modifications play a role in regulating FUS function. In this study, we discovered that FUS was acetylated at lysine-315/316 (K315/K316) and lysine-510 (K510) residues in two distinct domains. Located in the nuclear localization sequence, K510 acetylation disrupted the interaction between FUS and Transportin-1, resulting in the mislocalization of FUS in the cytoplasm and formation of stress granule-like inclusions. Located in the RNA recognition motif, K315/K316 acetylation reduced RNA binding to FUS and decreased the formation of cytoplasmic inclusions. Treatment with deacetylase inhibitors also significantly reduced the inclusion formation in cells expressing ALS mutation P525L. More interestingly, familial ALS patient fibroblasts showed higher levels of FUS K510 acetylation as compared with healthy controls. Lastly, CREB-binding protein/p300 acetylated FUS, whereas both sirtuins and histone deacetylases families of lysine deacetylases contributed to FUS deacetylation. These findings demonstrate that FUS acetylation regulates the RNA binding, subcellular localization and inclusion formation of FUS, implicating a potential role of acetylation in the pathophysiological process leading to FUS-mediated ALS/FTD.

ATAC-Me Captures Prolonged DNA Methylation of Dynamic Chromatin Accessibility Loci during Cell Fate Transitions.

DNA methylation of enhancers is dynamic, cell-type specific, and vital for cell fate progression. However, current models inadequately define its role within the hierarchy of gene regulation. Analysis of independent datasets shows an unanticipated overlap between DNA methylation and chromatin accessibility at enhancers of steady-state stem cells, suggesting that these two opposing features might exist concurrently. To define their temporal relationship, we developed ATAC-Me, which probes accessibility and methylation from single DNA library preparations. We identified waves of accessibility occurring rapidly across thousands of myeloid enhancers in a monocyte-to-macrophage cell fate model. Prolonged methylation states were observed at a majority of these sites, while transcription of nearby genes tracked closely with accessibility. ATAC-Me uncovers a significant disconnect between chromatin accessibility, DNA methylation status, and gene activity. This unexpected observation highlights the value of ATAC-Me in constructing precise molecular timelines for understanding the role of DNA methylation in gene regulation.

The Acetylation of Lysine-376 of G3BP1 Regulates RNA Binding and Stress Granule Dynamics.

Stress granules (SGs) are ribonucleoprotein aggregates that form in response to stress conditions. The regulation of SG dynamics is not fully understood. Permanent pathological SG-like structures were reported in neurodegenerative diseases such as amyotrophic lateral sclerosis. The Ras GTPase-activating protein-binding protein G3BP1 is a central regulator of SG dynamics. We found that the lysine 376 residue (K376) of G3BP1, which is in the RRM RNA binding domain, was acetylated. Consequently, G3BP1 RNA binding was impaired by K376 acetylation. In addition, the acetylation-mimicking mutation K376Q impaired the RNA-dependent interaction of G3BP1 with poly(A)-binding protein 1 (PABP1), but its RNA-independent interactions with caprin-1 and USP10 were little affected. The formation of G3BP1 SGs depended on G3BP1 RNA binding; thus, replacement of endogenous G3BP1 with the K376Q mutant or the RNA binding-deficient F380L/F382L mutant interfered with SG formation. Significant G3BP1 K376 acetylation was detected during SG resolution, and K376-acetylated G3BP1 was seen outside SGs. G3BP1 acetylation is regulated by histone deacetylase 6 (HDAC6) and CBP/p300. Our data suggest that the acetylation of G3BP1 facilitates the disassembly of SGs, offering a potential avenue to mitigate hyperactive stress responses under pathological conditions.

Deacetylase activity of histone deacetylase 3 is required for productive VDJ recombination and B-cell development.

Histone deacetylase 3 (HDAC3) is the catalytic component of NCoR/SMRT corepressor complexes that mediate the actions of transcription factors implicated in the regulation of B-cell development and function. We crossed Hdac3 conditional knockout mice with Mb1-Cre knockin animals to delete Hdac3 in early progenitor B cells. The spleens of Hdac3F/-Mb1-Cre+/- mice were virtually devoid of mature B cells, and B220+CD43+ B-cell progenitors accumulated within the bone marrow. Quantitative deep sequencing of the Ig heavy chain locus from B220+CD43+ populations identified a defect in VHDJH recombination with a severe reduction in productive rearrangements, which directly corresponded to the loss of pre-B cells from Hdac3Δ/- bone marrow. For Hdac3Δ/- B cells that did show productive VDJ rearrangement, there was significant skewing toward the incorporation of proximal VH gene segments and a corresponding reduction in distal VH gene segment use. Although transcriptional effects within these loci were modest, Hdac3Δ/- progenitor cells displayed global changes in chromatin structure that likely hindered effective distal V-DJ recombination. Reintroduction of wild-type Hdac3 restored normal B-cell development, whereas an Hdac3 point mutant lacking deacetylase activity failed to complement this defect. Thus, the deacetylase activity of Hdac3 is required for the generation of mature B cells.

ALS mutant SOD1 interacts with G3BP1 and affects stress granule dynamics.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease. Mutations in Cu/Zn superoxide dismutase (SOD1) are responsible for approximately 20 % of the familial ALS cases. ALS-causing SOD1 mutants display a gain-of-toxicity phenotype, but the nature of this toxicity is still not fully understood. The Ras GTPase-activating protein-binding protein G3BP1 plays a critical role in stress granule dynamics. Alterations in the dynamics of stress granules have been reported in several other forms of ALS unrelated to SOD1. To our surprise, the mutant G93A SOD1 transgenic mice exhibited pathological cytoplasmic inclusions that co-localized with G3BP1-positive granules in spinal cord motor neurons. The co-localization was also observed in fibroblast cells derived from familial ALS patient carrying SOD1 mutation L144F. Mutant SOD1, unlike wild-type SOD1, interacted with G3BP1 in an RNA-independent manner. Moreover, the interaction is specific for G3BP1 since mutant SOD1 showed little interaction with four other RNA-binding proteins implicated in ALS. The RNA-binding RRM domain of G3BP1 and two particular phenylalanine residues (F380 and F382) are critical for this interaction. Mutant SOD1 delayed the formation of G3BP1- and TIA1-positive stress granules in response to hyperosmolar shock and arsenite treatment in N2A cells. In summary, the aberrant mutant SOD1-G3BP1 interaction affects stress granule dynamics, suggesting a potential link between pathogenic SOD1 mutations and RNA metabolism alterations in ALS.

Effect of a Novel Engagement Strategy Using Twitter on Test Performance.

INTRODUCTION: Medical educators in recent years have been using social media for more penetrance to technologically-savvy learners. The utility of using Twitter for curriculum content delivery has not been studied. We sought to determine if participation in a social media-based educational supplement would improve student performance on a test of clinical images at the end of the semester. METHODS: 116 second-year medical students were enrolled in a lecture-based clinical medicine course, in which images of common clinical exam findings were presented. An additional, optional assessment was performed on Twitter. Each week, a clinical presentation and physical exam image (not covered in course lectures) were distributed via Twitter, and students were invited to guess the exam finding or diagnosis. After the completion of the course, students were asked to participate in a slideshow "quiz" with 24 clinical images, half from lecture and half from Twitter. RESULTS: We conducted a one-way analysis of variance to determine the effect Twitter participation had on total, Twitter-only, and lecture-only scores. Twitter participation data was collected from the end-of-course survey and was defined as submitting answers to the Twitter-only questions "all or most of the time", "about half of the time", and "little or none of the time." We found a significant difference in overall scores (p<0.001) and in Twitter-only scores (p<0.001). There was not enough evidence to conclude a significant difference in lecture-only scores (p=0.124). Students who submitted answers to Twitter "all or most of the time" or "about half the time" had significantly higher overall scores and Twitter-only scores (p<0.001 and p<0.001, respectively) than those students who only submitted answers "little or none of the time." CONCLUSION: While students retained less information from Twitter than from traditional classroom lecture, some retention was noted. Future research on social media in medical education would benefit from clear control and experimental groups in settings where quantitative use of social media could be measured. Ultimately, it is unlikely for social media to replace lecture in medical curriculum; however, there is a reasonable role for social media as an adjunct to traditional medical education.

HDAC6 regulates mutant SOD1 aggregation through two SMIR motifs and tubulin acetylation.

Histone deacetylase 6 (HDAC6) is a tubulin deacetylase that regulates protein aggregation and turnover. Mutations in Cu/Zn superoxide dismutase (SOD1) linked to familial amyotrophic lateral sclerosis (ALS) make the mutant protein prone to aggregation. However, the role of HDAC6 in mutant SOD1 aggregation and the ALS etiology is unclear. Here we report that HDAC6 knockdown increased mutant SOD1 aggregation in cultured cells. Different from its known role in mediating the degradation of poly-ubiquitinated proteins, HDAC6 selectively interacted with mutant SOD1 via two motifs similar to the SOD1 mutant interaction region (SMIR) that we identified previously in p62/sequestosome 1. Expression of the aggregation-prone mutant SOD1 increased α-tubulin acetylation, and the acetylation-mimicking K40Q α-tubulin mutant promoted mutant SOD1 aggregation. Our results suggest that ALS-linked mutant SOD1 can modulate HDAC6 activity and increase tubulin acetylation, which, in turn, facilitates the microtubule- and retrograde transport-dependent mutant SOD1 aggregation. HDAC6 impairment might be a common feature in various subtypes of ALS.